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Chiral phase properties of finite size hadronic systems are investigated within the Nambu-Jona- 
Lasinio model. Finite size effects are taken into account by making use of the multiple reflection 
expansion. We find that, for droplets with relatively small baryon numbers, chiral symmetry restora¬ 
tion is enhanced by the finite size effects. However the radius of the stable droplet does not change 
much, as compared to that without the multiple reflection expansion. 
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I. INTRODUCTION 

The behavior of (finite lumps of) quark matter is of great interests in cosmology, neutron stars, cosmic 
ray physics and heavy ion collisions Gill Absolutely stable nonstrange quark matter contradicts 
ordinary nuclei consisting of neutrons and protons. However, the existence of stable strange quark 
matter is still an open question and it may be realized in the form of strangelets, small lumps of strange 
quark matter. The (meta)stability of nonstrange and strange quark matter has been investigated within 
the MIT bag model 0, quark mass density dependent model 0J, Nambu-Jona-Lasinio (NJL) model 0, 
and so on j6j. In the MIT bag model, which assumes that asymptotically free quarks are confined in a 
bag, the bag constant and the current quark masses are phenomenological input parameters. Farhi and 
Jaffe 0 found a reasonable range of these parameters in which strange quark matter is stable, while 
nonstrange quark matter is unstable as compared to a gas of 56 Fe. The quark mass density dependent 
model takes account of the confinement mechanism by introducing density dependent quark masses. 
Fowler et al. |i] concluded that the critical density for a phase transition between nuclear and quark 
matter is sensitively dependent on the confinement mechanism. In the NJL model, the bag constant 
and constituent quark masses are generated dynamically in marked contrast to the above mentioned 
models. In Ref. [8f, it has been argued that stable nonstrange quark matter lies in chirally restored 
phase and, then, the model shows similar behavior to the MIT bag model (see, also, Ref. 0 ). The 
stability of strange quark matter has been also studied and concluded that strange quark matter is not 
absolutely stable when realistic values of current quark masses and coupling constants are used [l0]. 

On the other hand, unlike bulk systems, finite size effects cannot be ignored in the study of finite 
lumps of quark matter. So far, the MIT bag model and Fermi gas model including finite size effects 
have been widely used to study their properties 7, HH |H [12J, [lj, [Hj ■ For finite quark matter, it is 
known that the lumps are not energetically favored by finite size effects. The disfavor also arises from 
large strange quark masses. Therefore, finite size effects on dynamical quark masses (i.e. the behavior 
of chiral symmetry) might have significant effects on their properties. In this paper, we adopt the NJL 
model that is a simple tractable model respecting chiral symmetry. Then we consider isospin symmetric 
( n u = rid with rq = ( qjqi )) quark droplets. Finite size effects are taken account of by making use of the 
multiple reflection expansion (MRE) nn in in which is an approximation for the density of states in 
finite systems. The MRE itself has been derived for free massless particles or nonrelativistic particles. 
However, the MRE has been used to calculate the thermod yna mic quantities of massive quarks and 
reproduced well the results of the MIT bag model 0, UH Hj, HU • For instance, it can be checked that 
the vector and scalar (ifrp) densities in the MIT bag model are very similar to those computed 

in the MRE with little r dependence for bags of radius R > a few fm. Therefore, we expect that the 


* Electronic address: kiriyama@rcnp.osaka-u.ac.jp 
'Electronic address: hosaka@rcnp.osaka-u.ac.jp 



2 


MRE also works in the NJL model in the mean field approximation with a constant m inside the bag 
and apply the MRE to the effective potential of the model. 

This paper is organized as follows. In Sec. II, we formulate the effective potential of spherical quark 
droplets that consist of an equal number of up and down quarks. The Schwinger-Dyson equation and 
related thermodynamic quantities are also derived. In Sec. Ill, we present numerical results. Section 
IV is devoted to conclusions. 


II. THE MODEL 


In this section, the effective potential and related thermodynamic quantities of spherical quark droplets 
are derived in the two-flavor NJL model. We consider up and down quarks to be massless. Then, the 
SU(2 ) l x SU(2) fl chirally symmetric Lagrangian is given by 

C = + G [(4>ip) 2 + ('#Y5'?V>) 2 ] , (1) 

which is our starting point. Here, ip denotes a quark field with two flavors (Nf = 2) and three colors 
(N c = 3), and G is a dimensionful coupling constant (mass~ 2 ). The Pauli matrices t act in the flavor 
space. 

Let us first discuss the effective potential for bulk systems. In the mean-field (Hartree) approximation, 
the effective potential u> at finite temperature T and quark chemical potential /i is written as Q 


m 


f k 2 dk 


( 2 ) 


where v = 2 NfN c , Ek = \Jk 2 + m 2 , (3 = 1/T, with m being the dynamically generated (constituent) 
quark mass. Since the model is not renormalizable, we have to specify a regularization scheme. Through¬ 
out this paper, we use a sharp cutoff A = 600 MeV in the three dimensional momentum space with 
the coupling constant GA 2 = 2.45 sgj. With these parameters, the vacuum values of the pion decay 
constant / ff , the constituent quark mass m and the quark condensate are obtained as f n = 93.9 MeV, 
m = 400 MeV and (uu) = (dd) = (—244.9MeV) 3 . Further discussions about the stability of nonstrange 
quark matter can be found in Refs. 

We now turn to discussions on finite size droplets. Within the framework of the MRE in mm, 
we write the density of states for a spherical system as fc 2 pMRE/(27r 2 ), where 


67i 

Pmre = PMRE(fc, m,R) = 1 + -jTftfs 


12tt 2 
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( 3 ) 


with R being the radius of the sphere. The functions fs(k/m ) and fc(k/m) represent the surface and 
the curvature contributions to the fermionic density of states, respectively. The functional form of the 
surface contribution is given by 


fs 



1/2 k\ 
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( 4 ) 


On the other hand, the curvature contribution for an arbitrary quark mass has not been derived in the 
MRE. In this paper, we adopt the following Ansatz by Madsen fill : 


fc 
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Note that Eqs. 0 and © have the following to —> 0 limits: 

lim fs(k/m) = 0, lim fc(k/m ) = 

m —>0 m —>0 


( 5 ) 
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Thus, the effective potential for the spherical system which is regarded as a function of m , p and R is 
given by 


w 


m z 
4 G 
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(7) 


where the Fermi momentum, kF, is related to the chemical potential by p = yJkp + m 2 . Henceforth 
we restrict ourselves to zero temperature. For computation of a finite system, we choose a fixed baryon 
number A and a radius of the sphere R. The quark mass to is then determined by the Schwinger Dyson 
equation (SDE) that is the extremum condition of to with respect to to: duj/dm = 0. The resulting 
SDE is given by 


to = 2Gv 


d 


dm 


f A k2dk v 
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In addition to the SDE, since the baryon number is fixed to A , we have another equation, i.e. 

” r kF k 2 dk 


Vn B = V- 


2n 2 


'Pm re = A, 


( 8 ) 


(9) 


where V = 47tI? 3 /3 is the volume of the droplet and n B denotes the baryon number density in the 
droplet that is one third of the quark number density n q : 


n q If dui\ v f kp k 2 dk 

ns = y = “3^J T ^0 3 Jo yy^ MRE - 


( 10 ) 


Equations © and © have to be solved self-consistently. Therefore, by virtue of Eq. ©, one can reduce 
these equations to the following set of coupled equations: 


= 2 Gv 
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( 11 ) 

( 12 ) 


Notice that we consider kF in Eq. CD to be a variable independent of to. However, under the constraint 
CD kp depends on m. Therefore, these equations are also obtained by minimizing the energy density 


e 



k 2 dk 

2 ^_ 2 ~ PMHE-^k 


(13) 


with the constraint of fixed baryon number. 

Once we know the density dependent (i.e. the radius dependent) quark mass to, we can compute 
various physical quantities. For instance, the pressure p inside the droplet are given by 


P = 


TO 

4G 


I k F 


k 2 dk 3 A rz o 

-^-PMRE-C/fc + — Y fej, + TO 2 - p v ac, 


(14) 


where p vac has been introduced to ensure p —> 0 when kp —> 0 (i.e. R —> oo). 

We note that, unlike the bulk case, the SDE CD does not have a trivial solution to = 0. The reason 
is that pmre, the density of states of the MRE, depends on the quark mass to, and hence the r.h.s. of 
the SDE 1111) is no longer proportional to to, in contrast to the bulk case. 
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III. NUMERICAL RESULTS 

In this section, we present the numerical results. Unless stated otherwise, the baryon number A is 
fixed to A = 100. 

For a given baryon number, a nontrivial solution to the coupled equations CD and m is obtained 
as a function of the radius R. Then, as shown in Fig. 1, the pressure of chirally symmetric and broken 
phases are calculated from Eq. GJ. As a result, the critical radius at which the chiral phase transition 
occurs is identified as R), MRE) ~ 12 fm. The corresponding critical baryon number densit y n l MRE) is 
rather small n(, MRE) ~ 0.08no with no being normal nuclear matter density no = 0.17 fm - . [2fj 

The resultant quark mass as a function of the radius is shown in Fig. 2. Since the baryon number 
is fixed, by changing the radius R , the baryon number density ns changes. Therefore, Fig. 2 implies 
that the phase transition as ns is varied is of first-order when the effect of the MRE is included. For 
comparison, we also present the result without the MRE showing that the transition at R c ~ 4 fm is 
of second-order. The difference between these two curves obviously stems from the finite size effects. 
This result is roughly interpreted as follows. In the MRE, the finite size effects reduce the density of 
states and, therefore, the actual Fermi momentum increases for a given baryon number. The increasing 
tendency of the Fermi momentum is more pronounced when the baryon number is small (see Fig. 3). 
As a consequence, the chirally broken phase is energetically disfavored for small baryon numbers. 

Figure 4 shows the quark mass as a function of the radius for the cases of A = 100 and A = 1000. It 
is clear that, for relatively large baryon numbers, the first-order transition is weakened. In addition, at 
A = 1000 we find R(, MRE) ~ 14 fm and n(, MRE) ~ 0.5no- These results, off course, indicate that the finite 
size effects become less important for large baryon numbers. 

Let us consider the radius of the stable droplet from the usual pressure balance relation between the 
inside of the droplet and the vacuum (the outside of the droplet). The pressure inside the droplet with 
and without the MRE are shown in Fig. 5 for A = 100, where the vacuum pressure is taken to be 
zero. When R is small, the pressures are positive and chiral symmetry is restored inside the droplet. 
As R is increased the pressures decrease. Now without the finite size effects (MRE) the pressure keeps 
decreasing in the negative pressure region toward the cusp point, where the chiral phase transition takes 
place. Beyond the cusp point, the pressure curve follows the one in the chirally broken phase and turns 
into the positive region (but with small values), asymptotically approaching the zero value. There are 
two points where the pressure curve passes the zero, where the system is in equilibrium. The left point 
(R ~ 3.7 fm) corresponds to the absolute minimum of the energy per baryon, while the right point 
(R ~ 7.5 fm) corresponds to the local maximum. Consequently, the stable droplet of radius R s ~ 3.7 
fm lies in the chirally symmetric phase. 

On the other hand, with the MRE the pressure crosses zero only at R ~ 3.5 fm. We have calculated 
the pressure up to R = 10 8 fm and confirmed that the pressure of the chirally broken phase is always 
negative. Therefore, we find R< MRE > ~ 3.5 fm, which is somewhat small as compared to that without 
the MRE, and the stable droplet lies in the chirally symmetric phase. The radius of a stable droplet 
is presented in Fig. 6 as a function of A. As is evident from the figure, as long as A is sufficiently 
large, R ( S MRE) varies as A 1 / 3 . This result implies that, as mentioned before, the finite size effects are less 
important at large baryon numbers. Incidentally, we have observed the saturation property (R ~ A 1 / 3 ) 
when A is large. This is a consequence of the balance between the kinetic energy of quarks and the (bag 
like) volume type potential energy. 


IV. CONCLUSIONS 

We have studied the behavior of chiral symmetry in finite systems within the SU(2) L x SU(2) fl 
symmetric Nambu-Jona-Lasinio model. The multiple reflection expansion was used to take account 
of finite size effects. We considered spherical quark droplets with baryon numbers A = 100 — 1000. 
The Schwinger-Dyson equation for the quark mass was solved and, then, it turned out that the finite 
size effects enhance chiral symmetry restoration. The critical radii with and without the MRE differed 
significantly from each other. If the finite size effects are important also for strange quarks, it may have 
significant effects on the (meta)stability of strangelets. On one hand, the enhancement of restoring chiral 
symmetry lowers the constituent mass of strange quarks. This affects the argument made by Buballa 
and Oertel [h| who pointed out that the large strange quark mass at the relevant densities is the main 
reason that makes strange quark matter unstable. On the other hand, the finite size effects increase the 
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energy of strangelets BluLGHH. These two competing aspects should be carefully treated in realistic 
calculations 

In order to find the radius of the stable droplet, we calculated the pressure inside the droplet. Without 
the finite size effects (MRE) the pressure shows the characteristic behavior as shown in Fig. 5. Such 
behavior has been also obtained [3, |lj]] with interactions/parameters which differ from those in this 
paper. In contrast, with the MRE its behavior changes considerably. We found that the pressure of the 
chirally broken phase is always negative and, therefore, the chirally broken phase is unstable against 
collapse. However, the radius of the stable droplet does not change much as compared to that without 
the MRE. After all, it was found that chiral symmetry is restored in the stable droplet. This result 
reminiscent of the MIT bag model supports the former results without the MRE [gj, Q ■ 

We have already referred to the validity of the MRE in the first section. It should be also noted that 
the MRE contains several problems concerning its reliability. In the first place, pmre should have terms 
proportional to 1 /R 1 2 3 4 5 6 7 8 9 10 11 12 , 1 /i? , and so on, though their functional forms are not known. These terms would 
be dominant at small radii. Furthermore, it is well known that the MRE causes unphysical negative 
density of states at small radii. Neergaard and Madsen El have investigated the validity of the MRE 
and proposed a solution to these problems. However, for systems of A > 50, these problems have rather 
minor effects and the results presented here should hold when an improved density of states is employed. 

In the present analysis we have adopted the simplest NJL model with scalar and pseudoscalar channels. 
At the quantitative level the inclusion of vector and axial-vector channels could be important which 
brings repulsive correlations mmum. It might modify the present result such that, for instance, 
the transition point would be shifted to higher density side. However, we expect that our main conclusion 
drawn here on the relative effect of finite cavity as compared to the bulk case should be rather stable. 
This is, however, an interesting question to be investigated in the future work. 

Finally, we comment on the outlook for future studies. Physics of strangelets and the quark droplets 
at finite temperature should have important implications for the cosmological quark-hadron phase tran¬ 
sition and heavy ion collision experiments. It is worthwhile doing more detailed analysis including charge 
neutrality, chemical equilibrium by weak interactions and the color singlet constraint. Furthermore, it 
is also interesting to include other dynamical effects as well as the chiral phase transition, e.g., color 
superconductivity and color flavor locking j^j, [tI [24f. They may alter the present phase structure at 
small radii, because, at sufficiently low temperatures and high densities, the chirally broken phase would 
undergo a phase transition into the color superconducting phase. 
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FIG. 2: Radius dependence of the quark mass for A = 100 (solid line). For comparison, the result without the 
MRE is presented by dot-dashed line. 




FIG. 3: The density of states k 2 p^ ire/(2tt 2 ) (solid line) and that without the MRE (dotted line) in unit of 
MeV 2 as a function of k (left panel). The increase of the Fermi momentum due to the finite size effects Afci? as 
a function of A (right panel). The quark mass m and the radius R are taken to be their typical values; m = 200 
MeV and R = 10 fm. 
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FIG. 4: The radius dependence of the quark mass for the cases of A = 100 and A = 1000. For relatively large 
baryon number, the first-order transition is weakened. 




FIG. 5: The pressure inside a droplet (A = 100) as a function of R. The curves correspond to the pressure with 
(solid line) and without (dotted line) the MRE. Without the MRE the pressure becomes positive for R > 7.5 
fm, though it is difficult to see in this figure. 









